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ABSTRACT 

Aims. We study the expected variability patterns of blazars within the two-zone acceleration model putting special emphasis on flare 
shapes and spectral lags. 

Methods. We solve semi-analytically the kinetic equations which describe the particle evolution in the acceleration and radiation 
zone. We then perturb the solutions by introducing Lorentzian variations in its key parameters and examine the flaring behavior of the 
system. We apply the above to the X-ray observations of blazar lES 1218+304 which exhibited a hard lag behavior during a flaring 
episode and discuss possibilities of producing it within the context of our model. 

Results. The steady-state radio to X-rays emission of lES 1218+304 can be reproduced with parameters which lie well within the 
ones generally accepted from blazar modeling. Additionally, we find that the best way to explain its flaring behavior is by varying the 
rate of particles injected in the acceleration zone. 
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1. Introduction 

Blazars, a subclass of Active Galactic Nuclei, show variability 
which is observed across the electro-magnetic spectrum and can, 
in some cases, be as fast as a few minutes (Aharonian et al. 2007; 
Albert et al. 2007). The working blazar scenario is that we ob- 
serve their radiation coming from a jet which is directed, within a 
small angle, to our line-of-sight. The observed spectrum, which 
is clearly non-thermal in origin, shows a characteristic double- 
peak in a vFy plot, and is produced from a population of rela- 
tivistic particles, presumably accelerated by shock waves within 
the jet. While there are still open questions concerning the origin 
of the high energy gamma radiation, there is a consensus that the 
emission at lower frequencies, which usually extends from radio 
to X-rays, comes from electron synchrotron radiation. The ob- 
served variability can then be attributed directly to conditions in 
the acceleration region (for a recent review of the various blazar 
models see Bottcher 2010). 

The detailed X-ray observations of many blazars reveal that 
these objects show a rich and complex structure in their flar- 
ing behavior. Especially interesting are the trends in the lags 
between the soft and hard X-ray bands. The first observations 
which were detailed enough to show time structure in various 
energy bands (Takahashi et al. 1996) revealed the blazar Mrk 
421 having a soft lag flare - that is the soft X-rays lagged behind 
the hard ones. This was explained as synchrotron cooling: lower 
energy electrons cool with a slower rate than higher energy ones 
and therefore soft X-rays appear after the harder ones, i.e. soft 
lags hard. 

Subsequent observations, however, showed cases of the op- 
posite trend (Takahashi et al. 2000). This could not be explained 
by synchrotron cooling, but it was attributed to particle accel- 
eration. As fresh particles are accelerated to high energies, they 
radiate first soft and later hard photons, i.e. hard lags soft. 

Long, uninterrupted observations (Brinkmann et al. 2005) 
revealed that Mrk 421 shows both kinds of lags, with typical 
timescales of the order of lO^s. The same mixed behavior was 



also verified for the same source by Tramacere et al. (2009). 
Finally, blazar lES 1218+304 showed a hard lag during a giant 
flare (Sato et al. 2008). 

While several time-dependent models have been put for- 
ward to explain time variations (Mastichiadis & Kirk 1997; Li 
& Kusunose 2000; Bottcher & Chiang 2002), hard lag behav- 
ior needs a particle acceleration scheme to be taken explicitly 
into account (Kirk et al. 1998; Kusunose et al. 2000). Kirk et al. 
(1998, hereafter KRM) have shown that both soft and hard lags 
can be obtained during a flare once a detailed scheme for particle 
acceleration is considered. Adopting a two-zone model in which 
particles are energized in an acceleration zone by the first-order 
Fermi process and then escape and radiate in a radiation zone, 
they argued that flares produced by an increase in the number 
of the injected particles can produce soft lags if the radiating 
electrons are in the energy regime where cooling is faster than 
acceleration and hard lags when the two timescales are compa- 
rable. 



More recently, two-zone models which include SSC losses 
have been developed in Graff et al. (2008); Weidinger & Spanier 
(2010). In particular, Weidinger & Spanier include both first and 
second-order Fermi acceleration in their acceleration zone. A 
detailed study of the two types of acceleration and of their ef- 
fect on SSC spectra, but in an one-zone model, can be found in 
Katarzynski et al. (2006). 

In the present paper we investigate further the expected vari- 
ability within the two-zone model by letting its key parameters 
vary in time. While we keep the same assumptions as in KRM, 
we add also to the photon spectrum radiation coming from the 
acceleration zone. In Sect. 2 we present the basic points of the 
model and then examine some interesting cases in Sect. 3. In 
Sect. 4 we apply our results to lES 1218+304 and discuss them 
in Sect. 5. 



1 



K. Moraitis and A. Mastichiadis: X-ray variability patterns in blazars 



2. The model 

We employ the two-zone acceleration model as developed by 
KRM. According to this, electrons are accelerated by a non- 
relativistic shock wave which propagates along a cyUndrical jet 
and then they escape in a wider region behind the shock where 
they radiate the bulk of their energy. Following the KRM nota- 
tion we shall refer to the region around the advancing shock as 
the 'acceleration zone' (AZ) and to the escape region as the 'ra- 
diation zone' (RZ). The total size of the source is assumed to 
have a finite extent. 



2.1. Acceleration zone 

The evolution of the electron distribution function (EDF) in the 
AZ, Niy, t), is given by the continuity equation 



dN d N 

+ ^(Nj) + — = Q6(y - Tinj). 
at ay ?esc 



(1) 



The energy of each particle in this region changes with the rate 



y 2 

y=. ay 

tacc 



(2) 



where the first term describes first order Fermi-type accelera- 
tion with the energy independent timescale t^cc and the second 
describes the energy losses due to synchrotron radiation in the 
magnetic field B with 



a 



(3) 



Particles are injected in the AZ at some low energy ymjmeC^ with 
the rate Q and escape from it to the RZ with the energy indepen- 
dent rate t^^^. 

In our model we assume that variability can be produced by 
the change of one or more of the parameters face, fesc, Q and 
B in the AZ and thus treat them as time-dependent in Eq. (1). 
We define the dimensionless functions that include the time- 
dependency through the relations facc(0 = fa.o/aCO^ fesc(f) = 
fe,o/e(f), Q(t) = Qofqit) and B(t) = B()/B(f), where fa,o, fe.o, go 
and Bo are the unperturbed values of the parameters. 

The solution of Eq. (1) can be obtained semi-analytically 
with the method of characteristics. This method requires the so- 
lution of the initial value problem 



dy 
dt 



■ yiy,t), y(t*) = rinj. 



which can be written in the form 

y Tmax Tinj Tmax 

where 

Tmax = (ao4,o)~' 
and the functions ^a> are given by the relations 



vAt) = exp 



and 



Urtlt^ _dt_\ 
'Uhfini 





(4) 

(5) 
(6) 

(7) 
(8) 



The solution of Eq. (1) can then be written in the generic 
closed form 



N(y,t) = 



5(r;rinj,ri«) (9) 



where ^ = 1 -h fa,o/fe,Oj S{x;a,b) - \ if a < x < b and zero oth- 
erwise, and the function is given by Eq. (7) with /a replaced 
by /e. The quantity t, = t^{y, t) in the case where both the ac- 
celeration timescale and the magnetic field are time independent 
(/a = /b = 1) has the analytic expression 



tt = t — fa,0 If 



ymax/yinj 



Tmax/r ■ 



t) 



(10) 



while in every other case it is obtained from the solution of the 
impUcit Eq. (5). The upper limit of the EDF is obtained from 
Eq. (5) for u = and it is given by 



nW = 



(11) 



In the case where all parameters are time-independent, the 
unperturbed distribution function of the AZ is 



\s-2 



Noiy, t) = Qot^o {- —] y-' (l - — ) 

\ yinj ymax / \ ymax / 

5(r;rinj,ri(0) (12) 

and its upper limit is 



yiit) = 



1 + :^(e'/'"0 - 1)' 



(13) 



relation valid whenever /a = /b = 1 • The time it takes for steady- 
state to establish is a few tens fa,o and then the EDF is a power- 
law of index s up to the maximum Lorentz factor -/max (Eq. 6) 
where acceleration balances energy losses. 

2.2. Radiation zone 

The escaping particles enter the radiation zone where they ra- 
diate a part of their energy in the constant magnetic field Bq. 
Assuming that synchrotron losses dominate, we can write an 
equation for the evolution of their difi'erential density n{y, x, t) 
intheRZ 



dn d I 2 \ 



N{y,t) 



6{x - Xsh(O). 



(14) 



where Xsh(0 = Msh? is the position of the shock and Msh is the 
shock speed. The solution of Eq. (14) is easily found 



niy, X, f) = — — - — 5 (x; 0, Xsh(O) 

Mshresc(x/Msh) \ y 



with 



y. = 



(15) 



(16) 



Following KRM, we impose the additional restriction x > 
Xsh (f) - L on the spatial limits of Eq . ( 1 5 ) , in order to take into ac- 
count the finite size of the RZ, L. This is better expressed through 
the time the shock needs to travel the RZ 



th = L/Ufih- 



(17) 
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Particles not fulfilling this restriction are assumed to escape the 
RZ in a region of practically zero magnetic field and thus do 
not contribute to the total radiation. The EDF in the RZ is then 
obtained by spatially integrating the function n(y, x, f). 

In the case where all parameters are time-independent, the 
steady-state distribution in the RZ is a broken power-law hav- 
ing the same energy limits as the EDF in the AZ. The power- 
law has an index s for energies less than the breaking energy 
Till - ymaxfa,o/fb ^ud (s + 1) for y > ybr- The break is formed 
because of the finite size of the source, since the particles with 
y < 7br leave the source before they cool. Radiation from these 
particles is ignored because of the assumption that the magnetic 
field dechnes substantially outside the source. The time it needs 
the EDF in the RZ to reach steady-state is of the order of ?b- 

2.3. Radiation spectra 

Knowing the electron distribution functions in both zones, it is 
easy to compute the emitted synchrotron spectrum in the source 
frame through the relation 

= I dyN(y,t)Iy(y)+ \ dylyiy) \ dxn(y,x,t) (18) 
Ji Ji Jo 

where lyiy) is the single particle emissivity for synchrotron radi- 
ation that is given by 

V3e^B / V \ 

with 

Fix) = [Kii^{x)Kii^{x) - y [Kli^ix) - Kl^ix))j (20) 

(Crusius & Schhckeiser 1986) and the characteristic frequency 
is 

3 eB 

vo^- . (21) 

In Eq. (18) one can see the two differences that our approach 
has with the one adopted by KRM. The first is that we calculate 
radiation from both zones, while KRM ignore the contribution 
from the AZ. The second difl'erence is that we do not treat Ught- 
travel effects inside the source. This assumption is equivalent to 
taking Msh ^ c and in that case the shock speed is no longer a 
parameter of the problem. Our model thus considers the case of 
non-relativistic shocks. In the opposite case, light-travel efi'ects 
can be important depending on the frequency (see the discus- 
sion in KRM) and then the intrinsic variations are smoothed out 
on the light-crossing timescale. A study of these phenomena in 
the internal-shock model for blazars can be found in Bottcher & 
Dermer(2010). 

3. Flare profiles 

The model that we presented in the previous section can produce 
a wide variety of flaring behaviors by varying one or more of its 
basic parameters. The temporal behavior of these is essentially 
another free parameter as no theory can provide this information. 
While in many applications a short duration impulsive change 
is assumed (Mastichiadis & Kirk 1997, KRM), here we adopt 




13 14 15 16 17 18 19 

logv [Hz] 



Fig. 1. Snapshots of the synchrotron spectrum in the source frame. The 
curve labeled 'a' is the steady-state spectrum, while labels 'b' to 'g' are 
the snapshots at f = + 4fcfa.o for ^ = 2, 7 respectively. The vertical 
arrows define the frequency range (Vf , Vi) for which the light-curves of 
Fig. 2 are calculated. 

the Lorentz profile for the changes which has a pulse shape and 
analytic expression 

m to, w,n) = l + in- 1) ^ . (22) 

4(f - fo)^ + 

Ax t = to the pulse shows a maximum or a minimum, i.e. 
fi.{h;tQ,w,n) - n, depending on whether n > 1 or n < 1, 
while for \t - ?ol » ^72 it is /l = 1. The quantity w is the 
full width at half maximum (or minimum, FWHM) of the pulse, 
since /l(?o ± f Uo, w, n) = {n+ l)/2. 

3.1. Ciiange in tlie injection rate 

A change in the injection rate can be attributed to the encounter 
of the shock with a higher or lower density area. The increase of 
the injection rate leads to an increase in the overall normahzation 
of the EDFs and thus to the production of a flare. The case of a 
step-function change was examined in KRM. Here we study it 
for a Lorentz type change in more detail. 

A few snapshots of the synchrotron spectrum when the injec- 
tion rate varies according to fq{t) = /l(?; to, w, n) with w = 3?a,o> 
n = 3 and ?o » ?b - so that steady- state has been reached long 
before f = fo - are shown in Fig. 1. For the steady-state spec- 
trum, the set of parameters ymj = 10, B = 0.5 G, /b = 100/a,o> 
Qo = lO'^^s"! and ^a.o, h.o such tiiat ymax = 10* and s = 1.8, 
is adopted. The maximum synchrotron frequency in the source 
frame, v^ax = Tmax^o, is then a few keV. The corresponding 
light-curves are shown in Fig. 2. 

It is evident from Fig. 2 that the form of the flares is a func- 
tion of frequency. In order to quantify this, we define for each 
flare the following parameters. If y(t) = Iy(t)/Iy(0) is the form of 
the light-curve normalized to the steady-state luminosity, /v(0), 
then we define the time at which the flare peaks fpk and its peak 
value jpk = yitpk). We also define the flux doubling time during 
the rise of the flare, tr, and the flux halving time during the decay, 
fd, through the relations v(fpk - fr) = .y(fpk + td) = (1 +>'pk)/2. The 
ratio tr/ti is then a measure of the time symmetry of the flare 
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Fig. 2. Flares that are produced by a Lorentzian change of the injection 
rate in the frequency range vt = lO'^Vmax < v < lO^^Vmax = vi with step 
10''''vmax (from bottom to top, the first and the last curve are labeled 
with their corresponding frequency). 



and Wfl = fr + fd is the FWHM of the flare. We should note here 
that we do not fit each flare with a specific function, but rather 
we calculate the above quantities independently from the exact 
form of the flares. 

In Fig. 3 we plot the quantities fp^, fr/fd, vvg and >'pk as func- 
tions of frequency. In the first plot for v < 0.02vniax the time 
at which the flare peaks decreases with frequency, while for 
V > 0.02vn,ax it increases. This means that for v < 0.02vnjax 
the low-frequency flares precede the high-frequency ones, or, in 
other words, the flares show a 'soft-lag' . The opposite happens 
for V > 0.02vmax where the flares show a 'hard-lag'. 

As noted in KRM, who also find this result, soft-lags can 
be attributed to electrons for which the acceleration timescale is 
much faster than the cooling timescale, therefore for all prac- 
tical reasons these can be considered as pre-accelerated or in- 
jected at high energies with a ready power-law (see, for example 
Mastichiadis & Kirk 1997). On the other hand, hard-lags are ob- 
served at higher frequencies where the two timescales are com- 
parable, therefore the observer sees the wave of freshly acceler- 
ated particles moving to high energies. 

Further insight can be gained by considering the exact 
changes which the electron populations experience. The distri- 
bution function of the AZ can be written in this case as 



N{y, t) = Noiy, to, w, n) 



(23) 



with tt given by Eq. (10). This result indicates that the initial 
pulse propagates in the EDF from low to high energies, since its 
peak is located at y = yi(f - fo), as can be easily found by setting 
f, = fo- Moreover, this moving pulse preserves the shape of the 
initial pulse and thus it is symmetric around its center and has 
an amplitude equal to n. So, if we consider radiation only from 
the AZ, we expect to get roughly symmetric flares of the same 
amplitude which exhibit hard-lags. 

In the radiation zone the situation is more complicated, since 
the EDF there results from the integration of Eq. (15) along the 
extent of the source. The hard-lag behavior of the AZ is seen here 
also, but only in high energies, where the contribution from the 
particles residing just before the shock dominates in the integral. 
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Fig. 3. Characteristics of the flares shown in Fig. 2. 

In lower energies the perturbation propagates towards the oppo- 
site direction, as the intrinsic variations of N(y, t) are smoothed 
out by the spatial integration and only the effect of the fading 
pulse in the AZ stands out. 

In the second and third plot of Fig. 3 one sees that the width 
of the flares decreases with frequency, while the ratio fr/fj in- 
creases. This means that the flares become narrower and more 
symmetric as the frequency increases. If the rise time of the 
flares was connected with the acceleration timescale and the de- 
cay with the cooling timescale, then the ratio t^lti wiU be pro- 
portional to 



^acc 
^cool 



4,0 



r 



(aor) 



V 
^max 



1/2 



(24) 



The plot of this curve in the trlU diagram reveals that this is 
not the case and thus the relation of tr, ti to the acceleration and 
cooling timescales is more complicated. 

Finally, in the last plot of Fig. 3 one sees that the amplitude of 
the flares increases with frequency until it saturates to a constant 
value which is somewhat lower than the amplitude of the elec- 
tron variation. The level of saturation depends on the length of 
the pulse, because the longer the duration of the burst becomes, 
the closer the EDF gets to a steady state. 

3.2. Change in the escape timescale 

As a next case we examine a perturbation in the escape 
timescale. More precisely we assume that some physical pro- 
cess impedes momentarily particle escape, i.e. decreases the rate 
at which particles escape from the acceleration region. An in- 
spection of the EDF (Eq. 12) shows that an increase in the escape 
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Fig. 4. Flares that are produced by a Lorentzian change of the escape 
timescale in the frequency range Vf = 10"'vmax < v < 10''^Vmax = vi 
with step iC^'vniax (from bottom to top, the first and the last curve are 
labeled with their corresponding frequency). 



timescale leads to a decrease in the electron index. This will pro- 
duce a flare, the form of which, at various frequencies, is shown 
in Fig. 4. We assume that the change of the escape timescale fol- 
lows again a Lorentzian variation, i.e. /e(0 = fi^{t; to, w, n) with 
w - 3fa,o, « = 4/3 and to » t\„ while the parameters of the 
steady-state spectrum are the same as before. 

The specific features of the flares are shown in Fig. 5. As be- 
fore the flares exhibit a soft-lag in low frequencies and a hard-lag 
close to the maximum frequency. The frequency range in which 
hard-lags are observed is much smaller than in the previous case 
and also the flares now peak closer to the peak of the perturba- 
tion pulse, i.e. closer tot = to. 

The shape of the flares is asymmetric in time at low fre- 
quencies, having a rise time which is shorter than the decay one 
(second plot of Fig. 5). This however becomes symmetric, i.e. 
tt ^ td, at higher frequencies. The rising part of the curve indi- 
cates that tj/td is related to ?acc/^cooi and thus tiie rise of the flares 
is controlled by the acceleration and the decay by the cooling 
timescale. 

The width of the low-frequency flares decreases as in the 
previous case but in higher frequencies it increases. We note, 
however, that wg » w, i.e. even a narrow pulse in the escape 
timescale produces a rather broad X-ray flare. The amplitude of 
the flares increases with frequency as in the previous case, but 
now it does not saturate to a maximum value. 

The EDF of the AZ, Eq. (9), reads in fliis case 



N(y, t) = Noiy, t) exp 



w{n - \){s - 1) 



tan 



2V^ 
, t-tQ 
w yfn 



- tan 



w yjn 



(25) 



where /, is given by Eq. (10). The perturbation in the EDF prop- 
agates again from low to high energies but is no longer pulse- 
shaped and has a much broader extent. Moreover, the perturba- 
tion moves twice as fast now, since the location of its peak is 
given by y = 'yi(2(f - fo))- One notices that, qualitatively speak- 
ing, the trends in the present case are similar to the previous one. 
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Fig. 5. Characteristics of the flares shown in Fig. 4. 



There is a simple explanation for this. The perturbation intro- 
duced causes the escape timescale to increase first and then de- 
crease back to its original value. This in turn causes the electron 
spectral index to flatten and then steepen again. Therefore, the 
number of accelerated particles increases and then decreases to 
its unperturbed value. Thus this case can be considered as equiv- 
alent to the previous one which treated variations in the injection 
of particles at low energies. Note that, in the present case, the 
flattening of the EDF is partly compensated from the fact that, at 
the same time, the particles decrease in the RZ (see Eq. 15). 



3.3. Change in the acceleration timescale 

As a next example we investigate changes in the acceleration 
timescale and more specifically the variations induced when this 
decreases. This change has two effects on the EDFs and, there- 
fore, to the produced spectrum. First, an inspection of the EDF 
(Eq. 12) reveals that if the acceleration timescale decreases, then 
the electron index becomes flatter. Second, the maximum elec- 
tron energy will increase - see Eq. (6). The flares which are pro- 
duced for a Lorentzian change with w = 3/a,o and n = 3/4 are 
shown in Fig. 6. 

As frequency increases, the flare shapes become more asym- 
metric and peak at earlier times. In this case we get only soft-lag 
flares. This can easily be explained since the increase of ymax 
in the EDFs leads to the simultaneous increase of the maximum 
emitted frequency in the spectrum. This produces a large flare at 
high frequencies, since during pre-flare this part of the spectrum 
was in the exponentially decaying synchrotron regime. 
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Fig. 6. Flares that are produced by a Lorentzian change of the accelera- 



tion timescale in the frequency range Vf = 10 



< V < 10" 



vi with step 10 



(from bottom to top, the first and the last curve are 




W,0 

Fig. 7. Flares that are produced by a Lorentzian change of the magnetic 
field in the frequency range Vf = lO^'^Vmax < v < lO'^'^Vmax = V\ with 



step 10" 



(from bottom to top, the first and the last curve are labeled 



labeled with their corresponding frequency). 



with their corresponding frequency). 



frame is made through the relations 



3.4. Change in the magnetic field 

As a last case we examine variations of the magnetic field 
strength which are assumed to occur in the AZ only. A decrease 
in the magnetic field leads to an increase of the maximum elec- 
tron energy. Note that the maximum photon frequency increases 
also, even if B decreases and this results to a production of a 
flare. The form of the flares which are produced for a Lorentzian 
change of the magnetic field with w - 3fa,o and n - 0.8 is shown 
in Fig. 7. A new feature that appears is a small decrease of the 
flux preceding the low frequency flares. This is because the ini- 
tial drop of the magnetic field leads to the decrease of the total 
power radiated by the electrons and this is seen in low frequen- 
cies since the high-frequency spectrum rises due to the increase 
of Tmax which compensates for this reduction. 

The flares exhibit only soft-lags for a similar reason to that 
of the previous case. The flares also become narrower and more 
symmetric with frequency. In total, this case has many similari- 
ties with the change of the acceleration timescale, since in both 
cases the maximum electron energy follows the same temporal 
behavior. 



4. Application to 1ES 1218+304 

lES 1218-1-304 is a high-frequency peaked BL Lac object at a 
redshift z - 0.182. It was observed in May 2006 with Suzaku 
and the analysis of the observations (Sato et al. 2008) revealed 
that the flares in the various X-ray energy bands exhibit a hard- 
lag. Moreover, the flare shapes are asymmetric in time and their 
amplitude becomes larger with photon energy. These features 
resemble the ones found in the previous section for the cases of 
change of the injection rate or the escape timescale. 

The radio to X-rays spectrum of lES 1218-1-304 is shown in 
Fig. 8 together with the steady-state spectrum of our model as 
this is given in Eq. (18). The transformation to the observer's 



^ = <5^(1 + z) 



4ndl 



Vobs 



1+Z 



(26) 
(27) 



where 6 is the Doppler factor and = 880 Mpc is the lu- 
minosity distance of lES 1218-1-304, for the flat universe with 
cosmological parameters Oa = 0.7, Qm = 0.3 and Hq = 
70kms-' Mpc"^ 

The maximum photon frequency in the spectrum is ^ 
3.6 10^^ Hz, leading to the relation 



(28) 



We choose a random combination of these parameters so that 
Eq. (28) is satisfied and also take yinj = 10, fh = 20fa,o and 
s - 1.8. The value of fb is dictated by the frequency where the 
spectrum appears to break, v'^'^ ^ 9 lO'^ Hz, while the value of 
the electron index is chosen so that the photon index in radio fre- 
quencies is Ot ^ 0.6. The parameter yinj has a small eff'ect in the 
normalization of the spectrum and thus its value is unimportant. 

The remaining parameters can be further constrained by the 
variability timescale of lES 1218-H304. From the values of tp^ie) 
given in Sato et al. for the time at which the flares peak, we 
deduce that the variation responsible for producing these flares 
is the change in the injection rate. This is because the energy 
range in which /pk is rising is quite broad and rules out the case 
of the escape timescale change. As in Sato et al. we define the 
amount of hard-lag simply by the difference of the peak times of 
the flares in the various energy bands from the peak time of the 
flare in the highest energy band, namely 



Thard(e) - fpk(emax) - tpk(s). 



(29) 



where fimax = VSO keV is the logarithmic mean energy of 
the 5 - 10 keV energy band. The theoretical curve for Thard(£) 
for the Lorentzian change of the injection rate with parameters 
w = 1.7/a,o and n = 1.4, is shown in the upper panel of Fig. 9 
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Fig. 8. Radio to X-rays spectrum of lES 1218+304 (starred symbols, 
data taken from Sato et al. 2008) together with the steady-state theoret- 
ical spectrum (solid line) for the parameters given in text. The contri- 
butions from the acceleration zone (dashed line) and from the radiation 
zone (dotted line) are also shown. The optical data can be attributed to 
the host galaxy of lES 1218+304 (Riiger et al. 2010) and are thus not 
included in the fit. 

together with the observational data. From the fit we estimate 
the acceleration timescale 4,0 =1-9 10^ 610 s, where 610 = 5/10. 
This relation when combined with Eq. (28) enables us to ex- 
press the magnetic field and maximum electron energy of lES 
1218+304 as B = 0.06 6^^'^ G and = 1.3 10*' 6'^^^'\ 

Since the Doppler factor cannot be determined, we choose 
610 = 1 and the unperturbed value of the injection rate is then 
Qo = 3.2 10"'^ s"^ With these parameters we can estimate the 
size of the source L = 1.1 10'*'(Msh/0.1c) cm and the energy 
content of the source in electrons as £e = 3 lO'*^ erg. 

In the lower panel of Fig. 9 we plot the theoretical curve of 
the asymmetry of the flares, tr/ti, together with the observational 
data. Our definition of and ta differs from the one in Sato et al. 
but the ratio of the two quantities is the same in both cases. One 
sees that we reproduce the general trend of the data quite well but 
we overestimate the ratio ti/ta in the energy interval 1-3 keV. 

5. Conclusions 

In the present paper we have examined the expected variability 
of blazars in the context of the two-zone acceleration model. By 
varying four of the key parameters of the model, we found that in 
all cases flaring can be induced. In all cases the simulated flares 
exhibit soft lags - only in two cases, i.e. when varying flie injec- 
tion rate of particles or their escape rate, one can find frequency 
regimes where hard lags appear. Thus, we find that this scheme 
caimot produce hard lags at frequencies much lower than the 
maximum synchrotron frequency. Another conclusion is that the 
flares tend to mimic the electron variations at high frequencies. 
As a rule at low frequencies flares tend to be asymmetric with 
rise timescales much faster than decay ones. 

We were able also to make an acceptable fit to the hard lag 
flare of lES 1218+304 (Sato et al. 2008) by varying flie rate of 
particles injected into the acceleration process. We note that the 
fit can be improved if we were to vary two parameters instead of 
one but this lies outside of the aims of the present paper. The fit- 
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Fig. 9. Amount of hard-lag and the ratio tjtd of lES 1218+304 in the 
various X-ray energy bands together with the corresponding theoretical 
curves. 



ting we obtained gives parameters well within the ones accepted 
from blazar modeling. Note that the acceleration timescale de- 
duced (fa,o - 2 lO^^ios) might be long for first order Fermi accel- 
eration (Tamtni & Duffy 2009), however we have to emphasize 
that from the setup of the model, both acceleration and escape 
timescales were assumed to be independent of the energy y and 
the magnetic field strength B, therefore we do not try to connect 
them directly with theories of particle acceleration (Drury 1983; 
Blandford & Eichler 1987). 

We have restricted our analysis to the cases where syn- 
chrotron losses dominate and did not consider the inverse 
Compton scattering as a possible electron energy loss mecha- 
nism. This can be justified for sources where the energy den- 
sity due to magnetic fields exceeds the one due to photons. In 
this case the variabihty patterns due to inverse Compton scatter- 
ing will, in general, follow the ones due to synchrotron (Kirk 
& Mastichiadis 1999). On the other hand, an example of a 
model deaUng with X- and gamma-ray variability where in- 
verse Compton losses might be of importance can be found in 
Mastichiadis & Moraitis (2008) - note, however, that this is an 
one-zone acceleration model. 

Our model, despite its hmitations, is capable of producing 
a wide variety of flaring behaviors that could, in principle, be 
tested against the growing X-ray data on blazar variability. We 
note that recently, Garson et al. (2010) have also attributed the 
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properties of X-ray flares on intiinsic changes of the acceleration 
process. 
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